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We perform a model-independent phenomenological analysis of experimental data on proton form 
factor ratio. We find that only one of the two-photon exchange amplitudes, which we call SQm, 
is responsible for the discrepancy between Rosenbluth and polarization methods. The linearity 
of Rosenbluth plots implies that SQm is approximately linear function of e. The slope of SQm is 
extracted from experimental data and is shown to be consistent with theoretical calculations. 

At present, proton form factors (FFs) are measured with two techniques: Rosenbluth or longitudinal-transverse (LT) 
separation and polarization transfer (PT). The latter technique gives only FF ratio G e /Gm- When the experiments 
\ by the PT method were extended to Q 2 > 2 GeV 2 , the results have shown clear deviation from those of Rosenbluth 
method [l[. Two-photon exchange (TPE) effects are generally believed to be the cause for this discrepancy. They 
j were studied both phenomenologically @, [H, 0] and theoretically 0, Hi 0i IE ■ The theoretical calculation of TPE 
effects is nontrivial and technically difficult. On the other hand, assuming that the discrepancy is entirely due to TPE, 
one can extract some information on TPE directly from experiment. The values of TPE amplitudes, obtained in such 
a manner in Ref.Q, are in serious disagreement with the theoretical estimates. In the present paper we extract the 
TPE amplitudes from experiment more correctly and demonstrate that they agree with theoretical calculations. 

The following notation is used: particle momenta are defined by e(k) + p(p) — > e(k') + p(f>'); K = (k + k')/2, 
P = (p + p')/2, v = 4PK, transferred momentum is q = p' — p, Q 2 = —q 2 , the proton mass is M; the virtual photon 

Csl ■ polarization parameter is e — „i+q'^m2+q*) ■ 

| In general case the amplitude of the elastic electron-proton scattering can be written as 

M = -qT U '^ u ■ & \Fvf - W,7l^- + F 3 K„Yj^) U (1) 
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(the terms proportional to the electron mass are neglected). Invariant amplitudes (also called generalized FFs) Ft 
are scalar functions of two kinematic variables, say, Q 2 and e. For our purpose it is convenient to introduce linear 
Oh, combinations 

CD 

•fi; Ge = F 1 -tF 2 + ^F 3 , g M = F 1 + F 2 + e^F 3 , (2) 

^ ■ 
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. where t = . In general, the amplitudes are complex, but their imaginary parts give negligible contribution to the 
observables considered here. Thus later on, speaking of the amplitudes, we understand their real parts. 

With amplitudes ((2]), the reduced cross-section for unpolarized particles can be written like the Rosenbluth formula 

(± - £ ) ^2 f " + Q 2 \ 2 da 2 2 2 

The amplitudes can be decomposed as 

g E (Q 2 ,e) = G E (Q 2 ) + 6G ( p(Q 2 ,e) + 6G E (Q 2 ,e) + 0(a 2 ) (4) 

and similarly for Qm\ G e and Gm are proton electric and magnetic FFs and 8G E M + 5Q Et M are TPE corrections of 

order a. Here SG^ \ , denotes the part of the correction, calculated by Tsai [T3|. All infrared divergence is contained 
in it. The reduced cross-section then takes the form 

a R = eG\ + tG 2 m + 2sG E Sg E + 2tG m ^Qm + 2eG E SG ( P + 2tG m 5G ( P + 0(a 2 ). (5) 
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The terms containing SG EM are always subtracted from the cross-section by experimenters as a part of radiative 
corrections, so published cross-sections are, dropping terms of order a 2 , 



a R = eG E + tG 2 m + 2eGeSGe + 2tGmSGm- 



(6) 



These cross-sections are used as an input for the Rosenbluth separation. 
Now we will analyze Eq. ([6]) in the case of r > 1 . 

Since the TPE amplitudes are proportional to a k j^, it is natural to believe that their magnitude is about 1%. 
The key point is that Gm is enhanced with respect to Ge by about a factor of /z ~ 3 (proton magnetic moment), and 
even more according to PT data, so we have 



tG 2 m » eG% » 2eG E SG E , tG 2 m » 2tG m SGm » 2eG E 5G E - 



(7) 



Therefore the term 2eGe$Ge is much smaller than three other terms and can be safely neglected. Instead, the term 
2tG mSGm can be comparable with sG\ and thus strongly affect the results of Rosenbluth separation. We take for 
granted that Rosenbluth plots are linear (recent analyses 0, H| confirm this), i.e. op is linear function of e. This 
necessarily leads to the conclusion, that SGm is also approximately linear; we write it as 



5Gm{Q\s) = [a(Q 2 )+eb(Q 2 )]G M (Q 2 ), 



(8) 



where a and b are of order a. The coefficient a yields only a small contribution to the large e-independent term in 
<jp and does not change e-dependent term, thus we also may neglect it. After that 



gr — T G 2 M + e(G 2 E + 2rbG 2 M ) 
and the FF ratio squared obtained by the Rosenbluth method is actually 
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The FF ratio obtained in the PT experiments (with the inclusion of TPE) can be written as 

Ge 
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(9) 
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where Y 2l = jjtt-^-- The quantity e ^ + ^ never exceeds 0.17 and Y 2l should be about 1%, thus the bracket is very 
close to 1. Also, since the corrections 5Ge and SGm amount to about 1%, the ratio Ge/Gm is close to Ge/Gm and 
anyway matches it within the experimental errors, which are typically about 5%. Therefore we may say that PT 
experiments yield true FF ratio 
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Substituting (fT!?)) into (fTU|) . we obtain the TPE correction slope 



b - — (R^T - R 2 PT ). 



(13) 



We fitted R 2 PT from Ref.[Tl( and R\ T from Ref.(l2| with polynomials in Q 2 of the third order to minimize the x 2 



function 



X 



E 



R \ - p (Qi) 



2R t AR, 

where ARi are experimental errors. We have obtained, for Rosenbluth data 

Plt{Q 2 ) = m~ 2 (1.0736 - 0.1864Q 2 + 0.0358Q 4 H 
with x 2 = 24.2 for 28 d.o.f. and for PT data 



0.0007Q 6 ) 



(14) 



(15) 



P PT (Q 2 ) = A i" 2 (1.1184 - 0.3256Q 2 + 0.0323Q 4 - 0.0014Q 6 ) 



(16) 
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FIG. 1: (Color online) Experimental data and fits to FF ratio obtained by Rosenbluth method ([iRlt, solid curve and squares) 
and by PT method (/ii?pT, dashed curve and circles) and prediction for Rosenbluth method with positrons (dash-dotted curves, 
la range). 
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FIG. 2: (Color online) Extracted TPE correction slope b(Q 2 ). The dashed curves indicate estimated errors. 



with x 2 = 6-6 for 15 d.o.f. (everywhere Q is in units of GeV and \i = 2.793). The results are shown in Fig[TJ From 
this we derive the coefficient b: 

b = -0.0101/Q 2 + 0.0314 + 0.0008Q 2 + 0.0005Q 4 . (17) 

The dependence of b vs. Q 2 is displayed in FigO The coefficient a cannot be derived from the experimental data. We 
see that b varies from 0.01 to 0.05; if the unknown coefficient a is negative and equals about —6/2 then the relative 
size of TPE correction SQm/Gm = a + eb will never exceed 0.025, in consistency with the assumptions made at the 
very beginning. 

From the theoretical point of view, the TPE corrections can be split into elastic and inelastic parts. The elastic 
part was first calculated in Ref. [5j using simple monopole parameterization of FFs and then in a model-independent 
way in Refs. [g, The inelastic part is at present modeled by resonances Q or by parton models @. In Fig[3]we 
plot the elastic part of SGm, calculated as described in Ref.Q and straight line according to Eq.® with b given above 
and a chosen arbitrarily. A qualitative agreement is seen at all values of Q 2 . The gap between the curve and the line 
is small and should be possibly attributed to inelastic contribution. This results differ from those of Ref. Q , where 
authors, in fact, made three arbitrary assumptions about the TPE amplitudes (namely, in the notation of Ref.0], 
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that Y27, \Ge\, and \Gm\ are all independent of e), and obtained the quantity Y27 approximately 5 times larger than 
theoretical estimates 
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FIG. 3: (Color online) Comparison of extracted (dashed lines) and calculated (solid curves) values of TPE amplitude SQm/Gm- 

If the positrons are used instead of electrons, the TPE corrections change their sign. Thus we have for the Rosenbluth 
FF ratio squared, measured in positron-proton scattering 

^!t = #-2t6 (18) 

(FigHJ where dash-dotted lines indicate estimated la bounds for fiR^x)- The values of R 2 LT strikingly differ from 
R 2 LT and become even negative for Q 2 > 3 GeV 2 . Therefore the Rosenbluth measurements with positron beam at 
Q 2 > 2 GeV 2 would be an interesting test for TPE effects, and, on the other hand, a new way to measure true FF 
ratio, since the TPE corrections fully cancel in the average of R 2 LT and R 2 LT . Precision study of positron-proton 
cross-section is already proposed, for Q 2 < 1.6 GeV 2 [l3| and for 0.5 < Q 2 < 3 GeV 2 [Ill- 
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